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We study the evolution of orbital patterns in ferromagnetic layered ruthenates due to the compe- 
tition of Coulomb interactions, compressive c-axis and orthorhombic distortions in the presence of 
a polarizing orbital field coupled to the angular momentum. By means of the exact diagonalization 
on a 2x2 cluster and a cluster embedded analysis where inter-plaquette interaction is treated on 
mean field level, we determine the ground- state phase diagram. Specifically, we demonstrate that, 
via the activation of two or three of t2 g local orbital configurations, an external field applied along 
different symmetry directions can lead to inequivalent orbital correlated states. Starting from an 
antiferro-orbital pattern, for the easy axis case an orbital ordered phase is induced, having strong 
next nearest neighbors ferro-orbital correlations. Otherwise, a field applied along the hard axis leads 
a reduction of local orbital moment in a way to suppress the orbital order. 

PACS numbers: 74.70.Pq, 75.25.Dk, 71.10.-w 



I. INTRODUCTION 

Large Coulomb repulsion, orbital degeneracy and lat- 
tice distortions are widely accepted to be key parameters 
in determining the intriguing and fascinating phenomena 
of transition metal oxide materials^—. Colossal magne- 
toresistance (CMR) in perovskitic manganites is one of 
the most significant examples of the subtle competition 
between spin-orbital-charge degrees of freedom at the 
borderline between different types of ordered states^. 
Moreover, when electrons in degenerate d shells localize 
in Mott-Hubbard or charge transfer insulators, unusual 
collective behaviors may emerge as a consequence of the 
orbital frustration on the spin degree of freedom. In this 
context, quantum effects associated with the orbital de- 
gree of freedom are expected to be of great relevance both 
for etpi 7 - and ti g systems^. 

Recently, the physics of Ca3Ru207 compound has at- 
tracted a lot of interest due to a novel type of observed 
CMR 10- - that is likely to be driven by orbital rather 
than quantum spin effects. This compound belongs to 
the Ruddlesden-Popper series Ca„ + iRu„03„+i, n being 
the number of Ru-0 layers per unit cell, and thus ex- 
hibits a bilayered structure. The crystal structure is or- 
thorhombic with significant distortions in the ac plane 
due to the tilting of Ru06 octahedra around the b axis 16 . 
In zero field, Ca3Ru2 07 undergoes an antiferromagnetic 
(AFM) transition at K while remaining metallic 

and then a Mott-like transition at Tm/=48 K with a dra- 
matic reduction (up to a factor of 20) in the conductivity 
for T<Tmj i10 ' 11i13 ' 14i:l7 - This transition is accompanied 
by an abrupt shortening of the c axis lattice parameter 
below T MT 12.d£i. We notice that in the AFM state, neu- 
tron scattering measurements indicate that the magnetic 
moments align ferromagnetically within the double layer, 
and antiferromagnetically between the double layers (A- 
type AFM structure) along the c axi o 19 ' 20 . 
Concerning the mechanism behind the CMR in 



Ca 3 Ru207, the physics involved is likely to be fundamen- 
tally different from that governing all other magnetore- 
sistive materials. Experiments reveal that the electronic 
transport is so deeply interrelated to the spin-orbital 
correlations that, if one applies a magnetic field along 
the magnetization easy-axis a — , there occurs a first or- 
der metamagnetic transition to a spin-polarized state, 
without leading to a full suppression of the insulating 
state. Otherwise, the most favorable magnetoconducting 
state is driven by the field directed along the hard-axis 
b, where most likely an orbital disordered phase regime 
is obtained, without considerably affecting the magnetic 
structure. Indeed, in this regime the magnetization does 
not change significantly, thus suggesting that magnetic 
scattering is not the primary mechanism responsible for 
CMRU-i2.i2 i S2. These transitions are also accompanied 
by anisotropic structural changes and modifications of 
the coupling between the electrons and the lattice ^ 15 ' 18 . 
Notably, when the field is applied along the low temper- 
ature hard axis, there occurs a sharp structural change 
that ties the CMR to an increase of the c-axis lattice 
parameter — . The joining togheter of CMR, structural 
changes and metal-insulator transition is a clear indica- 
tion of the interplay between spin, charge and lattice de- 
grees of freedom — . 

Earlier theoretical studies on the phenomenology occur- 
ring in Ca3Ru2C>7 have been devoted to understand the 
nature of the antiferromagnetic metallic stat e 24 ' 25 and 
the character of the metal-insulator transition 26 . The 
crucial role of orbital degrees of freedom and structural 
distortions has been recently recognized in the conduct- 
ing behavior of Ca3Ru20?2£ while it has been exten- 
sively analyzed in the single layered compound Ca2Ru04, 
where it has gained a lot of extra interest in the context 
of orbital-selective Mott transitions^—. 
Having in mind the experimental scenario above men- 
tioned for Ca3Ru2C>7, the nature of the ground state 
within the layered ferromagnetic insulating state turns 
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out to be a relevant issue to be addressed. In this case, 
since the spin and charge degrees of freedom can be pro- 
jected out, a suitable microscopic description has to ex- 
plicitly include the degeneracy of ti g orbital sector, the 
compressive c-axis and the orthorhombic distortions of 
Ru06 octahedra. Then, as for other correlated insu- 
lating states, we base the low energy description on an 
effective orbital only model in the regime of strong on 
site Coulomb repulsion. Due to the many competing 
couplings involved in the problem, the analysis of the 
orbital correlations within the different ground states is 
performed via a systematic approach. Firstly, the given 
model Hamiltonian is solved by means of exact diag- 
onalization (ED) technique on a 2x2 plaquette. This 
analysis provides an unbiased insight on the character 
of the ground-state configurations as the microscopic in- 
teractions are treated without any approximation and on 
equal footing. Then, considering an embedded cluster ap- 
proach where the intra-plaquette Hamiltonian is solved 
via ED while the inter-plaquettes interaction is treated 
within a self-consistent Mean-Field (MF) approximation, 
we extract the phase diagram for the allowed broken sym- 
metry configurations. 

There are two main aspects behind the general complex- 
ity of the problem we faced; i) the Heisenberg-like Hamil- 
tonian for pseudospin-1 that we considered is marked by 
directional exchanges that in two-dimensions, due to the 
peculiar connectivity via the oxygen ligands of the tig 
centers, lead to antiferromagnetic correlations that are 
highly frustrated; ii) the tetragonal and the orthorhom- 
bic crystal field terms act like effective orbital fields that 
do not commute between each other. In this respect, 
even considering the case of a S=l/2 antiferromagnetic 
Heisenberg Hamiltonian, one can notice that the ap- 
plication of two non-commuting fields would produce 
a problem that cannot be solved exactly even in onc- 
dimension^r— We underline that our orbital case refers 
to a spin- 1 like system thus it intrinsically contains extra 
degrees of freedom whose interplay with the presence of 
two-non-commuting fields adds extra difficulties and dif- 
ferent scenarios as compared to the spin-I/2 situation. 
In particular, the octahedral distortions lift the pristine 
orbital degeneracy and set orbital patterns that privi- 
lege exchange along the easy axis a. This corresponds to 
lowering the local number of active degrees of freedom, 
that is reduced to the only two orbital configurations that 
can be effectively connected along the selected direction. 
Given this scenario, we show that the application of the 
field along the hard axis restores the local orbital fluc- 
tuations among the whole ti g subspace, by reintroducing 
the extra configuration that was quenched by the balance 
between orbital exchange and crystal field amplitudes. 
Special attention is devoted to analyzing the conse- 
quences of activating two or three of the tig local orbital 
configurations on the evolution of orbital order. One of 
the main findings is that a polarizing orbital field desta- 
bilizes antiferro-orbital ordered phases, both when apply- 
ing the field along the easy or the hard axis. However, we 



show that depending on the axial direction, we get two 
fundamentally different mechanisms: for a field applied 
along the easy axis a, the effect of the field is mainly to 
reorient the orbital moments inducing ferro-orbital cor- 
relations in the direction of the transverse applied field. 
Due to the orbital configurations involved, this situa- 
tion still has the features of an orbitally ordered state, 
where nearest-neighbor antiferro-orbital correlations are 
strongly suppressed. On the contrary, in the hard axis 
b case, the activated configuration along b reduces the 
amplitude of the ordering orbital moments in zero field 
up to reach an orbital disordered phase. 
We finally consider possible connections between the 
symmetry breaking associated with the orbital ordered 
state and the peculiar magneto-conducting behavior of 
layered ruthenates. 

The paper is organized in the following way: in the next 
section we introduce and write down the microscopic 
model; in Sec. Ill we present, for a 2 x 2 plaquette, 
the numerical results obtained for the phase diagram as 
function of the microscopic parameters stressing the role 
played by the purely structural distortions and their in- 
terplay with external magnetic field while in Sec. IV we 
present embedded mean field calculations based on exact 
diagonalization of 2x2 plaquette; the last section contains 
a discussion of the results presented and the conclusions. 



II. THE EFFECTIVE ORBITAL MODEL 

Based on the assumption that the planar spin configu- 
ration is not substantially modified by the applied field, 
we consider the coupling between the field and the lo- 
cal orbital angular momentum contribution as the domi- 
nant one. In particular, starting from different orbital or- 
dered (00) ground states at zero field, we will study how 
these configurations are modified as a consequence of the 
induced orbital "flip" via the coupling between the ap- 
plied field and the orbital angular momentum along dif- 
ferent symmetry in-plane directions. Then, by means of 
these results, we will consider the evolution of the ground 
state under the external field, and in turn use such in- 
vestigation to discuss the possible implications with re- 
spect to the anomalous CMR phenomenon observed in 
the Ca3Ru207. 

As mentioned in the Introduction, to investigate the 
correlations developing within the orbital sector, we start 
from a situation where the system is fully polarized and 
all the spins are aligned along a given crystallographic di- 
rection, in accordance to the observed A-type configura- 
tion. In this frame, the attention is focused on the orbital 
degree of freedom, out of the competition between the 
Coulomb repulsion (orbital exchange), octahedral distor- 
tions, and the effect of the orbital field. In particular, 
having in mind the Ca-based ruthenates, the case of flat 
octahedral configuration with orthorhombic distortions 
is considered. In such a scheme we can explore the inter- 
play between orbital exchange and the quenching of or- 
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bital correlations provided by structural distortions. The 
specific model Hamiltonian we refer to is built up by the 
following contributions: 

H = H exc + H te t + H ort + H a , (1) 

where H exc stands for the orbital exchange interaction, 
H te t and H ort are the crystalline fields terms, describ- 
ing local c-axis tetragonal distortions and in-plane or- 
thorhombic deformations, respectively, and H a repre- 
sents the coupling of the orbital angular momentum to 
the external field. 

The first term in Eq. ([I} is derived from the superex- 
change Hamiltonian for S — 1 spins and it is able to 
describe, in the low energy limit, the interaction aris- 
ing from the virtual excitation between two neighboring 
sites. In a previous work 3 ^, we have proved that such 
a description, which has the same roots as the superex- 
change in a Mott insulator with nondegenerate orbitals 7 , 
can be applied to the case of ruthenates, being consistent 
with a significant value of an on-site intraorbital Coulomb 
interaction. 

In ruthenates, the Hund's coupling is larger than the 
crystal field splitting of ti g subspace and the Ru 4+ ions 
are in triplet configuration. Dealing with an atomic t\ g 
configuration, the superexchange interaction arises from 
virtual processes of the type dfdlj — > dfd 3 —> dfdlj. 
Each of the orbitals involved is orthogonal to one cubic 
axis; for instance, d x is orthogonal to the x axis while 
d y and d z are orthogonal to y and z, respectively. Due 
to the connectivity via intermediate oxygen ligands, only 
two out of the three ti g orbitals can be effectively con- 
nected for each axial direction. 

The spin/orbital model for Ca-ruthenates can be di- 
rectly derived from the one adopted for cubic vanadates 8 , 
by performing a particle-hole transformation within the 
tig sector, in a way to map the problem of four elec- 
trons to that of two holes and consequently the doublon 
states in 'no hole' configurations. In this way, we re- 
duce ourselves to consider virtual excitations of the kind 

d 2 t d] -> did] -> d?4 

Here, we refer to a simplified version of the Hamiltonian 
in Ref. 35 , obtained under the assumption Jh/U — > 0, Jjj 
being the Hund's coupling and U the ( large) Coulomb in- 
teraction. Moreover, in the present case, the spin degrees 
of freedom have been projected out, by assuming a fully 
polarized planar configuration, where (Si-Sj) = 1/4. The 
orbital exchange Hamiltonian is then written as follows: 

H exc = ±JJ2 E & -Ti + \n inj ]W , (2) 
t (y)ll7 

where f7 — {rj i5 r yi , rj^} are pseudospin operators, act- 
ing in the subspace defined by the two orbitals which can 
be linked along a given 7 direction. 

For our purposes, it is convenient to characterize the cor- 
relations of each orbital flavor in terms of the 'no hole' 
states x, y, z, standing for the local orbital configurations 
having no hole in the corresponding orbital d x , d y , d z . 



For instance, if the bond is along the z-axis, the virtual 
hopping connects x and y only, and the orbital interac- 
tion may be expressed via the Schwinger representation: 
T % 1 = P^iyPix, i±i = P^ixPiy, t c zi = i(n iy - rife), and 
nf = 2n iz + m y + n ix , where n, 7 = p^ 7 Pi 7 and and 
Pi-y are creation and annihilation operators for 7 configu- 
ration at site i. Similar relations with exchange of bond 
index hold for the other axis directions. Note that, to 
avoid any confusion with the pseudospin components, we 
renamed the orbital flavor z = c. Moreover, J = 4t 2 /U, 
t being the hopping amplitude assumed equal for all the 
orbitals in the ti g manifold, due to the symmetry rela- 
tions of the connections via oxygen it ligand o 35 ' 36 . 
It is worth to notice that, in absence of crystal field terms 
lifting orbital degeneracy, orbital exchange described by 
Hamiltonian in Eq. [2]is intrinsically frustrated: due to or- 
bital singlets composition, antiferro-orbital correlations 
along x preclude exchange along y and viceversa. The 
result is a superposition of valence bond states which are 
uncorrelated along orthogonal directions. 

Concerning the crystalline field contributions, H te t is 
the tetragonal term linked to the compressive or tensile 
character of the c-axis octahedral distortions. In the 'no 
hole' state representation, H te t is written as follows: 

Htet = A tet y^\n iz - -{n iy + n lx )] . (3) 

i 

Eq. ([3|) implies that negative A tet values lower the energy 
of z configurations, thus simulating a compressive distor- 
tion. This case is relevant for describing the short-c-axis 
insulating phase of Ca 3 Ru 2 07 

The H or t term describes in-plane orthorhombic deforma- 
tions as due to a difference between the a and b crys- 
tallographic axis. From a general point of view, the or- 
thorhombic deformations are introduced by means of the 
coupling between the orbital degree of freedom and the 
strain fields s\ = s xx — s yy corresponding to a rotation 
axis like [100] and £2 = s xy corresponding to a rotation 
axis [110] in the tetragonal basis [xyz]. These distor- 
tions tend to lift the orbital degeneracy and thus can be 
viewed as effective orbital polarizing fields. From symme- 
try considerations^ 3 ^, one can show that a general form 
for the orthorhombic term can be expressed as follows: 

H ort = X)(Ai j0 t£ + A 2 , r^) . (4) 

i 

where Ai. Q and A2 iQ correspond with the strain fields e\ 
and £2- Having in mind the orthorhombic distortions for 
the Ca3Ru2C>7 compound, where the RuOg octahedra are 
severely tilted in the ac plane thus rotated with respect 
to the b axis, one can assume that the A2.0 term is the 
relevant contribution to be included. Therefore, in our 
analysis the last term in Eq.(l) is just the following: 

H ort = A J2r c xl . (5) 
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Still, negative values of A D in H ort , tend to stabilize a 
quantum state given by a superposition ~ x + y. This 
corresponds to select a preferential axial direction for or- 
bital connectivity in the (a, 6, c) frame that is obtained 
by performing a (7r/4, 7r/4, 0) rotation of (x,y,z) in real 
space. In the following we will refer to this as the 'or- 
thorhombic' frame and to the corresponding orbitals as 
d a , dt,, d c . 

Being interested in the response to an orbital only cou- 
pled field, we also include in the Hamiltonian of Eq. (TTJ) 
the coupling between the external field and the local an- 
gular moment^, i.e. 

i/ Q = 2B^r«, (6) 

i 

where B is the magnetic field coupled to the orbital 
degrees of freedom only in unit of Bohr magneton and 
a — a, b defines the direction, along which B is applied. 

III. GROUND-STATE ORBITAL PATTERNS 
FOR A 2 x 2 PLAQUETTE 

In this section we present the ground-state (GS) di- 
agrams, obtained by means of an exact diagonalization 
study of the Hamiltonian above introduced on a 2 x 2 pla- 
quette. As a strategy, we follow the evolution of the GS 
induced by the c axis compressive distortions, orthorhom- 
bic deformations and the external orbital only field. The 
orthorhombic field is assumed in a way that the a direc- 
tion is the easy-axis for the angular moments. 

It turns out to be easier to characterize the different 
ground-state behaviors by evaluating the density and or- 
bital correlators in the orthorhombic frame. We therefore 
introduce 'no hole' configurations a, b, c, where the empty 
state is localized at d a , db, d c , respectively. In Fig.[TJb), 
graphical representation of a, b, c, together with original 
x, y, z (Fig. [TJ a) is reported. In terms of this new no- 
tation, it is easy to define the action of the pseudospin 
operators in the orthorhombic basis for the a, b, c flavors. 
For example, along the c-axis, the pseudospin in the or- 
thorhombic frame acts on b and a and it can be expressed 
via T% i = p\ b p la , t c _ 1 = p^ ia p ib , r c zi = \{n ih - n ia ). 

Finally, in order to simplify the reading of the phase di- 
agrams, we introduce the following notation: FO, AFO, 
WO stand for ferro-orbital ('no-hole' state on neighbors 
homologue orbitals), antiferro-orbital ('no- hole' state 
on neighbors off-diagonal orbitals) and weak antiferro- 
orbital correlations, respectively. Besides, P-FO stands 
for configurations with partial FO character, C-AFO for 
canted AFO while PO indicates the para-orbital high- 
field configuration. Furthermore, the superscripts (a, 6, 
c) stand for the direction of the main contributions for the 
orbital correlators in terms of the pseudospins operators. 
The underscript (x, y, z) refers to the main character of 
the orbital correlations in the pseudospin space, i.e. it 
provides indication of the anisotropy of the orbital pat- 
tern. A schematic view of the representative orbital pat- 



tern configurations that contribute to the ground state 
on the 2x2 plaquette is reported in Fig. [TJ c)-m) , where 
the main component of the ground state is represented 
by means of the local density distribution of 'no-hole' 
configurations, in a way that a specific orbital on a given 
site stands for the local configuration having no holes in 
that orbital. 

We notice that, except for the FO/AFO regions that are 
expected without explicit calculations, all the remain- 
ing phases are untrivial since they show an intermediate 
character, where the correlations have competing behav- 
iors. For completeness, here we specify the way those 
regions are classified: i) PO stands for the para-orbital 
uncorrelated high-field-configuration, driven by the field 
H Q , where the pseudospins are x a oriented and the cor- 
relators (t"t"), evaluated in the orthorhombic frame, 
have the maximum saturated value of 0.25 [see for ex- 
ample Fig. 10 a) and 12 b)]; ii) generic C-AFO stand 
for regions where nearest-neighbor antiferro-orbital cor- 
relations characterizing AFO are strongly suppressed and 
coexist with FO components, as shown for example in 
Fig 10. iii) WO has the feature of an isotropic weak 
antiferro-orbital phase, obtained when the field is applied 
along the hard axis b, resulting from a reduction of the 
nearest-neighbors (t^t^ ) antiferro-type correlations due 
to a non- vanishing onsite momentum t*. In this region, 
FO and AFO correlations are the same order of magni- 
tude, as shown in Fig. 12, meaning that no preferred 
behavior can be identified. 



A. Effect of structural distortions 

As a first step we discuss the case at zero applied field. 
In Fig. [21 the ground-state diagram is shown, obtained 
by varying the crystal field microscopic parameters A tet 
and A„ in unit of J. It turns out to be convenient to 
span the different regions of Fig. [5] by fixing an ideal 
horizontal line corresponding to e.g. |A Q /J| = 0.8, while 
varying the tetragonal parameter. 

Tetragonal crystal field acts as an effective orbital bias 
that tunes the density of c configurations from to 1, 
corresponding to weak and strong regime of compressed 
octahedra, respectively. Looking at the relevant cor- 
relation functions, we see that the two ending regions 
are both characterized by a complete quenching of the 
orbital correlations. In the weak tetragonal regime, the 
orthorhombic field is strong enough to give rise to a 
complete FO a (see Fig. [TJ panel c)), where stabilizing 
b configurations generates full FO correlations for the 
t" components. In the extreme flattened RuOg regime, 
instead, we get a FO c region (see Fig.[TJ panel c)), where 
the ferro-orbital correlations are saturated along the z 
direction for the flavor c. 

These two regions are linked by intermediate regions 
characterized by a partial filling of the the ti g subspace. 
We highlight that, due to the constraint of the planar 
geometry, and that of the orthorhombic field favoring b 
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FIG. 1. (Color online) Schematic representation of the possi- 
ble orbital patterns on a 2x2 plaquette, generated by the in- 
terplay between orbital exchange, tetragonal and orthorhom- 
bic distortions, applied orbital field. First panel shows the 
axis notation, panel a) and panel b) the graphical represen- 
tation of x, y, z and a, b, c configurations, respectively. Panels 
c)-m) report the planar orbital pattern in terms of the local 
distribution of 'no-hole' configurations; FO a /AFO a stand for 
orbital patterns having homologue configurations on neigh- 
bors sites, for a — a, b, c flavor. P-FO stands for configura- 
tions with partial FO character, C-AFO for canted AFO. PO a 
indicates 'para'-orbital states where where all pseudospins are 
Ty oriented; WO stands for configurations characterized by a 
softening of AFO a correlations, due to nonvanishing (t^) . 



configurations, the orbital interaction along 6, involving 
a and c, is inactive thus yielding an anisotropic effect 
between a and b directions. 

Particularly, in the regime of intermediate tetragonal 
field, 0.2< \A tet /J\ < 0.4, we find an AFO a z region, rep- 
resented in Fig. [T]e), characterized by the alternation of 
c and b, that gives rise to predominant AFO correlations 
along z direction, for flavor a. This region is interesting 
when compared with the magnetic/orbital pattern in 
Ca 3 Ru 2 07; indeed, the occurrence of an AFO pattern 
can explain, in a spin/orbital picture, the existence of 
ferromagnetic in-plane correlations, and it can be also 
related to the poorly metallic behavior of the GS, as due 
to an 00 pattern. 

We finally remark that, in the AFO", the order of mag- 
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FIG. 2. Diagram summarizing the behavior of the GS config- 
urations of a 2 x 2 cluster as a function of A /J and Amt/J- 
Stars stand for smooth crossovers in the relevant correlators, 
solid dots for abrupt jumps. Phases are labelled in accordance 
to Fig. [I] 



nitude of the x, y, z correlators for the a flavor is compa- 
rable, though the amplitude for the z channel is larger 
than the others. 



B. Interplay between structural distortions and 
orbital polarizing field 

1. Easy axis 

Let us now discuss the interplay between the structural 
distortions and the effects of a polarizing field. In doing 
that, we consider just one representative value for the 
orthorhombic strain |A Q /J| = 0.4. This is the most in- 
teresting regime where the structure of the ground state 
is more sensitive to external small perturbations. This 
would permit to explore the modifications induced on the 
corresponding line of Fig. [21 and particularly of the region 
at the boundary between AFO a and FO c , occurring for 
intermediate/strong tetragonal crystal field amplitude. 

Due to the anisotropy introduced by the orthorhom- 
bic field, we distinguish between the response along the 
a and b directions. The corresponding ground-state dia- 
grams are displayed in Figs. [3] and 01 The investigation of 
these diagrams will be confined to the low fields regime 
since, approaching high fields, the orbital momentum be- 
comes fully polarized, as in the PO a region of Fig. Q]h), 
resulting from stabilizing ~ (b + ic), or the PO h region of 
Fig. [1]1), corresponding to ~ (c — ia) local configurations. 
We do believe that, in this limit, the assumption of ne- 
glecting the spin magnetization becomes weak and thus 
an analysis of the complete spin-orbital effective model 
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FIG. 3. Ground-state diagram of a 2 x 2 cluster obtained for 
an orbital-only applied field along the easy axis a as a function 
of \A te t/J\, for a fixed value |A o /J|=0.4. Dots correspond to 
abrupt jumps in the evolution of the relevant off-site orbital 
correlators; stars and diamonds to smooth crossover in the 
correlation functions. Different regions are labelled according 
to Fig. [1] 



FIG. 4. Ground-state diagram of a 2 x 2 cluster obtained for 
an orbital-only applied field along the hard axis b as a function 
of \Atet/J\, for a fixed value |A o /J|=0.4. Dots correspond to 
abrupt jumps in the evolution of the relevant off-site corre- 
lators; stars and diamonds to smooth crossover. The regions 
are labelled according to Fig. [T] 



is required. 

Fig. [3] shows the occurrence of crossovers in the or- 
bital correlators among different GS configurations gen- 
erated for B\\a. We do notice a peculiar evolution of 
the boundary for the AFO" state. Indeed, aligning the 
moment along y for the a flavor, regions AFO a and FO c 
move towards a new regime where AFO a along y is par- 
tially reduced evolving into a canted like state C-AFO" 
(see Fig. [TJ panel i)), with a non-zero total angular mo- 
mentum, directed along y, coexisting with partial AFO a 
correlations along z. 



2. Hard axis 

The response for B\\b is instead more stiff (see Fig.[3|. 
In Fig.[5j we give a schematic representation of the action 
of pseudospin operators on different orbital subspaces: 
the ground state at zero field is a combination of b and 
c states and the application of the field along b induces 
a-type configurations that are unfavored with respect to 
the orthorhombic deformations. This implies that a much 
higher value of the field is required to activate the orbital 
exchange along b. Along this symmetry direction, the or- 
bital polarizing field tends to restore a description where 
all the three t2 9 flavors are active. Due to the new config- 
urations a, and to the competing mechanisms involved, 
the resulting effect is a weakly ordered WO configuration. 
Such state is marked by a coexistence between AFO cor- 
relations along a and FO along b (see Fig. [TJ panel m)). 



Now, comparing the response of the system for the differ- 
ent field directions in the interval of crystal field given by 
about 0.1< A tet /J < 0.5, one can observe that, if the 
field is along the a axis, the closest configurations that 
can be activated, starting from AFO a and FO c , still have 
an orbital ordered pattern (OO) as given by C/AFO a . 
On the contrary, for B\\b, due to the activation of the b 
flavor, we get a changeover from an OO pattern to an 
WO configuration. 

This analysis shows that when AFO and FO configura- 
tions are close in energy, the application of the external 
?field may act as a tunable parameter within the orbital 
configurations, thus suggesting that the complex field re- 
sponse in the bilayered ruthenate may be consequence of 



Ground State Ella (Easy axis) 



IS h (Hard axis) 

id 
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FIG. 5. (Color online) Schematic representation of orbital 
configurations activated by an applied orbital field. Left panel 
shows 6 and c orbital configurations set by the interplay be- 
tween crystal field amplitudes. Middle panel shows the orbital 
subspace where r a is active; right panel shows the orbital sub- 
space where r h is active. The orbital field, applied along hard 
axis b, restores orbital fluctuations among the whole orbital 
subspace by re-introducing a configurations. 
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FIG. 6. Schematic representation of a 2x2 cluster embedded 
in a periodic environment. Filled square denotes 2x2 pla- 
quette described by Eq. [T] whose interaction with the field 
produced by four surrounding empty squares is expressed by 
Eq. [7] Dashed square defines the basic unity of the periodic 
structure. 



the subtle interplay between different microscopic mecha- 
nisms that involve the orbital couplings and the fields as- 
sociated with the octahedral deformations in presence of 
large Coulomb coupling. One of the main findings is the 
occurrence of an orbital ordered pattern at the boundary 
between FO and AFO states, stabilized when the field 
is along the easy axis, as due to the effective freezing 
of one orbital flavor. When the field is applied along 
the hard axis, quantum fluctuations associated with the 
three-orbital degrees of freedom are restored and tend to 
destroy the orbital ordering. 



change interaction between the plaquette and the near- 
est neighbors sites in the {x, y) plane, treated on a mean 
field level. Here, the average (...) of the local pseudospin 
operators is determined after solving the full quantum 
problem within the plaquette. The decoupling of the or- 
bital exchange on nearest-neighbour plaquettes tends to 
prefer broken symmetry states with AFO correlations. 
This tendency is in competition with the local distortive 
fields due to the octahedral distortions and is strongly re- 
lated to the intra-plaquette quantum effects. Hence, the 
solution is determined in a self-consistent fashion, accord- 
ing to the following scheme: we assign initial conditions 
for (t^) and (nj), and use them to evaluate improved 
expectation values; the procedure runs until convergence 
is achieved, with the requested accuracy. This scheme is 
applied for all the different initial conditions associated 
with possible orbital states emerging from ED on the 4- 
site plaquette. In particular, we consider all the starting 
configurations where local pseudospins are oriented along 
a given direction and FO, A-FO, C-AFO and WO cor- 
related, and we finally choose the lowest energy orbital 
configuration among the ones self-consistently obtained. 
Mean field drives the solution towards broken symme- 
try states having a specific orbital pattern that can be 
characterized by local order parameters (f^), and by ex- 
plicitly evaluating orbital correlations developing along 
different axial direction of the orthorhombic frame. We 
point out that this feature is not intrinsic in the MF ap- 
proach but is strongly related to the anisotropic micro- 
scopic environment due to crystal field terms. The latter 
break rotational invariance in the pseudospin space, al- 
lowing solutions with a specific orbital patterns to be 
stabilized. 

We finally remark that, by exploring broken symmetry 
states with defined 00 patterns, we can get informations 
about the thermodynamic stability of the corresponding 
00 phases. 



IV. CLUSTER EMBEDDED CALCULATION 

In this section we extend the cluster calculation consid- 
ering the case of a 4-site plaquette embedded in a periodic 
environment as schematically represented in Fig. El Here, 
the intra-plaquette Hamiltonian is solved via exact diag- 
onalization and the inter-plaquettes interaction is treated 
within a self-consistent Mean-Field approximation. 

Indeed, the MF scheme of analysis corresponds to de- 
couple the inter-plaquette terms of the Hamiltonian as it 
follows: 

= j E E E WW) + i^i>] (7) - 

J E E E \{<){^) + \{n i ){n j )p\7) 

-y=x,y p,=x,y,z (ij)\\~t 

where 7 stands for the axial direction and fi for the pseu- 
dospin component. Hence, Eq. ([7]) represents the ex- 



A. Effect of structural distortions 

We first consider the MF evolution starting from or- 
bital configurations of a 2 x 2 plaquette at zero applied 
field. We look for most stable solutions, chosen among 
those having FO-AFO and C/AFO character. In Fig. |7J 
we show the phase diagram as a function of the tetragonal 
and orthorhombic crystal field amplitudes. When com- 
paring this diagram to the cluster calculation (FigfSJ) , the 
following features emerge: 

i) the value of \A tet /J\ separating AFO" and C-AFO a re- 
gions is not modified by the MF and scales with |A Q / J| 
according to the same linear law; 

ii) in the moderately/highly flattened region, the most 
stable solution has a C-AFO a character. Looking at or- 
bital correlations in the orthorhombic frame, this corre- 
sponds to AFO correlations in x, y components of flavor 
a coexisting with FO correlations for r°. Approaching 
the extremely flattened side, the AFO correlations for 
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FIG. 7. Diagram summarizing the behavior of the GS con- 
figurations as a function of A /J and At e t/J for a cluster 
embedded in a periodic environment. 



r° are gradually replaced by the FO character of r° r~, 
which is saturated reaching FO c . Moreover t£ is decreas- 
ing and FO-correlated. 

iii) In the regime of strong orthorhombic field, we observe 
that the full FO a is no more stable and leaves the place to 
a region having a uniform P-FO a character (the evolution 
of the expectation values as a function of the orthorhom- 
bic field suggests that FO a is stabilized for higher values 
of |A Q | > 2). This configuration is somehow symmetric 
with respect to C-AFO a , but has an inverted balance be- 
tween \z) and (x,y) sector. 

iv) in the AFO" the self-consistence leads to a state with 
pure z components of the a and b flavor (correlations 
along x and y renormalize down to zero while they were 
not vanishing for the uncoupled cluster case.) 



B. Effect of orbital-only applied field 

1. Easy axis 

The effect of a polarizing field along easy axis a is to 
induce a nonvanishing y component of local orbital mo- 
mentum (Ty). This anisotropic effect is thus relevant for 
the evolution of the C-AFO a , where AFO a correlations 
exist both for the x and y components. Therefore, we 
expect that the polarizing field naturally wil select the 
latter. Moreover it's interesting to see how the induced 
FO correlations compete with the general tendency to- 
wards AFO. In Fig- El we report the evolution of the line 
| A Q / J | = 0.4 of Fig. [3 as a function of the applied field. 
As one can see, one of the main differences, with respect 
to the 4-site plaquette calculation (see Fig- EJ) , is that now 
most of the transitions occur as a smooth crossover in the 
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FIG. 8. (Color online) Schematic representation of possi- 
ble orbital patterns for a 2 x 2 cluster embedded in a peri- 
odic environment as due to the interplay between orbital ex- 
change, tetragonal and orthorhombic distortions, and applied 
orbital field. Panels a)-f) report the planar orbital pattern 
in terms of the local distribution of 'no-hole' configurations; 
FO a /AFO a stand for orbital patterns having homologue con- 
figurations on neighbors sites, for a = a,b,c flavor. C-AFO 
for canted AFO. PO Q indicates 'para'-orbital states where 
where all pseudospins are Ty oriented; WO stands for config- 
urations characterized by a softening of AFO a correlations, 
due to nonvanishing (r^) 



relevant correlators, represented by empty dots and stars 
(see Fig. ED). 

High field regions are uninteresting (PO phases are sta- 
bilized), so we skip them in the discussion. In the low 
\A tet /J\ region, we see that the AFO" region remains 
quite stable and the effect of the polarizing field is to 
progressively orient the local orbital momentum from z 
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FIG. 9. Diagram summarizing the behavior of the GS con- 
figurations as a function of Atet/J and B a /J, for a cluster 
embedded in a periodic environment at |A /J| = 0.4. 
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FIG. 10. (Color online) Expectation value of orbital correla- 
tors in the orthorhombic frame as a function of applied field, 
evaluated for \Atet/J\ = 0.6 and |A /J| = 0.4. Panels a)-b) 
show dominant contributions arising when the field is applied 
along easy axis. 



to y. The main modifications occur in the moderately 
flattened region where the C-AFO region-that now has 
broken symmetry, being active only for y component 
leaves place to Ci-AFO a j,. This region has the features 
of a strengthened OO, where c and a configurations alter- 
nate along the bonds [see Fig. [5]d)]. Looking at orbital 
correlators in Fig.[TU]a) and b), we deduce that in this re- 
gion the competition between crystal field amplitude and 
applied field has the effect to suppress nearest neighbors 
exchange. The OO that is stabilized has a strong FO 
character along opposite diagonals for Ty and t°, respec- 
tively, as shown by next nearest neighbors correlations, 
which are almost saturated to a value of ~ 0.25. 
To clarify the mechanism leading to the suppression of 
AFO correlations, we point out that the characteristic 
energy scale is of the order of J. The analogy with 
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FIG. 11. Diagram summarizing the behavior of the GS con- 
figurations as a function of Atet/J and Bb/J, for a cluster 
embedded in a periodic environment at |A Q /J| = 0.4. 



an Heisenberg antiferromagnet, disordered by an ap- 
plied magnetic field, suggests that orbital correlations are 
tuned by the interplay between applied field and orbital 
exchange, in a way that is not dependent on the crys- 
tal field amplitude, once the orthorhombic deformations 
have selected the 'easy flavor' a. 



2. Hard axis 

When the field is applied along the hard axis b, a 
stronger intensity than that for the a case is required 
to stabilize the full PO region. The evolution towards 
PO is the result of the smooth suppression of AFO cor- 
relations, as demonstrated by the suppression of (Ty/Ty ) 
between nearest neighbors [see Fig.[T2la)]. In the regime 
0.4 < \A tet /J\ < 0.6, at the boundary between AFOJ 
and C-AFO", this determines a weakly antiferro-orbital 
(WO) behavior, represented in Fig. |8]e). This region is 
the result of disordering C-AFO a by introducing a non- 
vanishing onsite momentum (r*) along the hard axis di- 
rection. As we already pointed out in previous section, 
this corresponds to restore orbital fluctuations along b 
by introducing a configurations, as illustrated in Fig. [5] 
Reactivating fluctuations in the whole orbital subspace, 
brings the system back to a situation where frustration 
of the orbital exchange is enhanced. 
We notice that, contrary to the easy axis case, for the 
hard axis orbital exchange between nearest neighbors is 
not hindered. Moreover, we remark that the energy scale 
stabilizing WO is now strongly related to the crystal field 
amplitudes. We explored the evolution of region bound- 
aries in Fig. [TTJ and found that the onset of WO region 
grows linearly with A D . These aspects underline that the 



10 




0.2 0.4 0.6 0.8 1.0 




FIG. 12. (Color online) Expectation value of orbital correla- 
tors in the orthorhombic frame as a function of applied field, 
evaluated for \Atet/J\ = 0.6 and |A /J| = 0.4. Panels a)-b) 
show dominant contributions arising when the field is applied 
along hard axis. 



softening of AFO correlations is here caused by the reduc- 
tion of onsite orbital momentum due to the new activated 
flavor, instead of by the field-induced FO correlations. 



V. CONCLUSIONS 

We have here presented a study of field tunability of 
orbital correlations for layered ruthenates in an insulat- 
ing ferromagnetic configuration. We have analyzed the 
evolution of the possible orbital patterns under the effect 
of the external magnetic field coupled to the local an- 
gular momentum only. Numerical calculations on a 2x2 
cluster and of a cluster embedded in a periodic environ- 



ment reveal that the competing effects driven by orbital 
exchange and crystal field put the system on the verge 
of a FO/AFO behavior, a circumstance that makes the 
character of OO in absence of an applied field quite elu- 
sive. This outcome may be supported by recent Resonant 
X-Ray scattering measurements on CaaR^Oy— , where 
no appreciable signal of staggered FO or AFO is found 
within experimental resolution. We argue that this fea- 
ture is at the origin of the soft behavior and the easy 
field- tunability of the OO phases. Particularly, AFO ex- 
change leaves the place to strongnext-nearest neighbors 
FO correlations when an external orbital field is applied 
along the easy axis, while we get a tendency to a isotropic 
weak antiferro-orbital correlated state for the hard axis. 
An interpretation of these results can be summarized as 
follows: the orthorhombic deformations tend to remove 
the orbital degeneracy and lead to an ordered state where 
only two local t^g configurations are active. The effect of 
the polarizing field is to keep the orbital correlations in 
the two-orbital sector when the field is along easy-axis; 
otherwise, along the hard-axis the field forces the system 
to allow for all possible local orbital fluctuations. The 
latter situation is responsible for demolishing the OO 
ground state set by the competition between Coulomb 
and crystal field interactions. Speaking on a general 
ground, in many multi-orbital Mott insulators the metal- 
insulator transition occurs as a symmetry breaking phase 
transition leading to magnetic, orbital and charge order- 
ings. We speculate that a possible connection may exist 
between the demolition of long range OO by an applied 
field and a likely metallic behavior and that this connec- 
tion may be of interest when referred to unconventional 
magneto-conducting properties of Ca3Ru207. 
We finally point out some important remarks about a 
quantitative connection of our calculation to the phys- 
ical parameters regime expected in the layered ruthen- 
ates, and specifically of Ca3Ru207. According to our re- 
sults, the evolution of region boundaries between AFO, 
C-AFO or WO behaviors is mainly concentrated within 
the windows 0.2 < A tet /J < 0.8 and 0.1 < A Q /J < 0.4. 
Checking the range of the microscopic parameters that 
bias the proposed orbital patterns is not easy since the 
amplitude of J, A tet and A Q is difficult to be directly 
measured and require a model derivation. An explicit 
ab — initio evaluation for Ca3Ru20y is not yet avail- 
able. For the closest members of Ruddlesden-Popper 
series according to their physical and structural prop- 
erties, Ca2Ru04 and Si^R^Oy respectively, one ends 
up with the following ab — initio estimates J ~ 50meV, 
A te t ~ 20meV and A Q ~ 5meV— <Q. These values set 
the ratios A tet /J ~ 0.4 and A D /J ~ 0.1, whose orders 
of magnitude are compatible with our choices of analysis 
for the regions of phase diagrams where major changes 
occur. 
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